Holographic quantum states 
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We show how continuous matrix product states of quantum field theories can be described in terms 
of the dissipative non-equilibrium dynamics of a lower-dimensional auxiliary boundary field theory. We 
demonstrate that the spatial correlation functions of the bulk field can be brought into one-to-one corre- 
spondence with the temporal statistics of the quantum jumps of the boundary field. This equivalence: 
(1) illustrates an intimate connection between the theory of continuous quantum measurement and quan- 
tum field theory; (2) gives an explicit construction of the boundary field theory allowing the extension 
of real-space renormalization group methods to arbitrary dimensional quantum field theories without the 
introduction of a lattice parameter; and (3) yields a novel interpretation of recent cavity QED experiments 
in terms of quantum field theory, and hence paves the way toward observing genuine quantum phase 
transitions in such zero-dimensional driven quantum systems. 



In recent years we have witnessed tremendous success 
in the calculation of physical properties of quantum many- 
body systems from their wavef unctions. This development 
has been spurred by studies of the entanglement properties 
of strongly correlated quantum spin systems: it has been 
established that natural states of quantum lattice systems 
are only slightly entangled, and hence typically obey an 
entropy area law IBBl. These new developments have also 
shown that the manifold of physical quantum lattice states 
are well captured by matrix product states (MPS) or finitely 
correlated states (FCS) f5\. 

These developments have also allowed the interpretation 
of the renormalization methods of Wilson IS) and White 
Q |8l as applications of the variational principle to ma- 
trix product states, and have led to natural generalizations 
of these RG schemes to higher dimensions IIQ UTTI . A key 
property, crucially responsible for the success of these RG 
schemes, has been the fact that the information concern- 
ing the quantum correlations of a natural quantum lattice 
state is encoded in the variational parameters of an auxil- 
iary zero-dimensional system. 

A natural next step is then to develop a similar approach 
for quantum field theories, and this is exactly the subject 
of this paper. However, capturing the manifold of low- 
energy wavefunctionals in, e.g., bosonic theories is much 
more challenging due to the continuous infinity of degrees 
of freedom (as opposed to the lattice of finite-dimensional 
Hilbert spaces for quantum spin systems). The most natu- 
ral way to proceed is to discretize the continuous degrees 
of freedom by imposing a lattice cut-off and truncating the 
local hilbert spaces at each site |12-14|. This is similar 
to the approach taken by lattice gauge theory which, so 
far, provides essentially the only systematic way to under- 
stand nonperturbative effects in quantum field theories. Re- 
cently, however, it was established that there is no need to 
impose a lattice cutoff because continuum limits of ma- 
trix product states can be directly defined, and these states, 
termed continuous matrix product states (cMPS), can rep- 
resent the low-energy physics of non-relativistic field theo- 



ries extremely accurately llT5l[T6l . 

In this work we describe a method to generate cMPS 
for quantum field theories. Our procedure is based on the 
paradigm of continuous measurement ifTTl . where one sub- 
jects a quantum system to a weak sequence of measure- 
ments of some physical observable. This procedure is a 
natural generalisation of the sequential preparation scheme 
for MPS proposed in (2T\ to the continuous setting and we 
show that this procedure naturally generates cMPS of Ref. 
lITSl and we relate our approach to that of ||T2| . The per- 
spective offered here also allows one to design more flex- 
ible classes of states which can be tailored to a system of 
interest, including quantum states for bosonic systems at 
arbitrary filling. Our approach also suggests a natural gen- 
eralisation to arbitrary dimensions: in the case of a 2+1 
dimensional quantum field theory the resulting boundary 
field evolves according to local dissipative dynamics in l-i-l 
dimensions. An alternative interpretation of the resulting 
wavefunctionals is that the boundary field provides a local 
parameterisation of the bulk field; this realises one of the 
major prerequisites identified by Feynman for the success- 
ful application of the variational principle to quantum field 
theory |18|. 

Surprisingly, the bulk and boundary fields have a direct 
interpretation in the context of cavity electrodynamics for 
trapped atoms llT9lJ22l : the role of the auxiliary system is 
played by the trapped atom and the quantum field describes 
the photons leaking from the cavity. An atom with a fairly 
low number of internal addressable levels (e.g. D = 6), 
would already allow the reproduction of all static correla- 
tions functions in, e.g., the Lieb-Liniger model |23|. This 
is achieved by observing the temporal counting statistics of 
the photons leaking from the cavity. The present paper also 
sheds new light on the recently discovered phase transitions 
of the quantum trajectories obtained in dissipative systems 
ll24l : such dynamical phase transitions are in correspon- 
dence with static quantum phase transitions of a quantum 
field theory in one dimension higher. Just as classical one- 
dimensional nonequilibrium systems exhibit phase tran- 



sitions similar to their two-dimensional static analogues, 
0+1 dimensional nonequilibrium systems exhibit dynam- 
ical phase transitions analogous to quantum phase transi- 
tions in 1+1 dimension. 

We begin by modelling a measurement of some physi- 
cal observable M on a quantum system with dimension D, 
which we initially call the "system". Our model, known 
as von Neumann's prescription 0251 . is defined as follows. 
We attach a quantum system with a continuous degree of 
freedom, called the meter, in a fiducial state vector |0) and 
couple it the system for some time t according to the inter- 
action Hj = M ® p. Supposing the system is initially in 
I i?i>) , we see that after the interaction the state of the system 
and meter is given by 
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where \mj), j = 1,2, ... ,D, are the eigenstates of M 
with corresponding eigenvalues rrij, \(j)) = J2i=i </'jl"^j)' 
\x = rrijt) = D{a,^/2mjt)\x = 0), and D{a,a) = 
e"" ~"° is the (phase-space) displacement operator. In our 
case, at t = 1 the displacement operator simply effects a 
translation of the meter's state, initially localised at posi- 
tion x = 0, to the locations of the eigenvalues x = nrij. It 
is important to note that here and in the sequel we never ac- 
tually perform a projective measurement of the meter. This 
can be thought of as corresponding to the situation where a 
projective measurement of the meter is performed but the 
measurement record is discarded. While this is an admit- 
tedly rudimentary model of the physical measurement pro- 
cess it does afford a considerable potential for generalisa- 
tion. 

Now, the core of our proposal is to turn von Neumann's 
measurement prescription on its head and regard the meter 
(attached as an ancillary system) as the fundamental sys- 
tem A and the system as an auxiliary ancilla B. In this 
way we can think of it as a state generation device: we 
can obtain a variety of physical quantum states of the me- 
ter A alone — a bosonic system with a continuous degree of 
freedom — ^by exploiting the measurement prescription and 
then tracing out, or perhaps measuring, the system B. In 
these terms we have a way to generate states of a quantum 
system with a continuous degree of freedom. The chal- 
lenge remains, however, to somehow exploit this procedure 
to obtain quantum states of a continuous infinity of such 
continuous degrees of freedom. 

The way we do this here is to model the continuous 
measurement of a specific POVM (positive operator valued 
measure ll26l ). Our model, which follows Ref. 1 17| closely, 
is defined by a family of D x D complex matrices R{x), 
X G [0, L] IITtII which we instantaneously and infinitely 
weakly measure on B at time t = x, which is additionally 
evolving according to some free hamiltonian K{t). We do 
this by introducing a collection ^ of n meters, labelled by 



r = 1,2, . . . ,n. The total hamiltonian is given by 

H{t) = Kit)®I_^ + Hi{t), (1) 



where 
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We are interested in the limit where n — )• oo and e — )• 
with ne = L fixed. (Indeed, for finite n, this approach in- 
cludes the scheme of Ref. lfT2ll when locally representing 
the meters in terms of coherent states). The choice of the 
coefficient ^/e in the definition of Hj is motivated by gen- 
eral considerations; any other scaling would lead to trivial 
dynamics, thanks to the quantum Zeno effort, or to the sit- 
uation where the meters and system do not interact ifTTIl . 

We now supply an interpretation of the cMPS of lITSl 
based on the sequential preparation prescription of ll22l 
and a continuous measurement scenario ||28]| . It is 
straightforward to integrate the Schrodinger equation for 
Eq. |l|: U{L) = Te'^^o "^"^'^^ where T denotes time- 
ordering, so that U{L) = V{L, L - e)W{L - e)V{L - 
e,L- 2e)W{L - 2e) • • • W{e)V{e, 0), where V{b, a) = 

By making the standard definition of the discretised quan- 
tum field operator as ^(re) = arej \f(- we can, in the limit 
e — )• 0, describe this evolution via 

In our continuum limit the collection A of the meters may 
be regarded as a geometrically one-dimensional bosonic 
quantum field. The auxiliary system B may be regarded as 
a geometrically zero-dimensional bosonic quantum field. 

We now learn that the evolution Eq. (J2]) prepares cMPS. 
Indeed, if we initialise the meters A in the vacuum state 
vector |ri)_4 of the quantum field we can, for each r, exploit 
the Baker-Hausdorff formula to first order to arrive at the 
identity 

valid to 0(e), to rewrite the limit Eq. (J2]l for U{L){\Q)^<^ 

Ib) as 

where Q{s) = —iK{s) — ^B){s)R{s). This is — up to a 
trace over the auxiliary degree of freedom — identical to the 
definition of a cMPS [15|. 

Thus, a cMPS is a quantum state of the quantum field A 
preparable via a continuum measurement scenario: i.e., we 
initialise the quantum field A in some (known) prespeci- 
fied quantum state uja = |^).a(^| and adjoin an auxiliary 
zero-dimensional quantum field B initialised in some (pos- 
sibly mixed) fiducial state p. We then interact A and B ac- 
cording to the manifestly unitary continuous measurement 



dynamics U{L). We then discard the auxiliary system by 
tracing it out to obtain the quantum state 



(^A 
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We pause here to emphasise an important point: there is 
no reason for the initial state uja of the quantum field A to 
be the vacuum. All we require is that initial field state arises 
from the continuum limit of a product state cOr 'S> ■ ■ ■ 
LOi of the meters r = 1, . . . , n. Thus, crucially, we can 
allow for initial field states with a high density of bosons 
and superpositions of bosons. This will be important in 
a variety of contexts, particularly those pertaining to dense 
systems with nonlinear interactions, where a cMPS defined 
using the vacuum will be insufficient. Physically this can 
be explained as follows: the interaction Hj between A and 
B can only transport e bosons to ^ in a time e. Thus the 
interaction can never achieve superpositions of terms with 
more than a constant number of bosons per unit length. 

Suppose o"^ is a cMPS. We now investigate the dynam- 
ics of the auxiliary system throughout the continuous mea- 
surement process by instead tracing out the field A. Again, 
we first consider the discrete setting and take the contin- 
uum limit. We set p(0) = ps and 

p{re) = tv_A [U{re)iujA ^ p{0))W{re)] , 

where now U{re) = Te^^Ja ^(*)''*. We then consider 
-{p{{r + l)e) — p{re)) and expand U{re) to second or- 
der (just as it is done when describing dissipative quantum 
systems in the weak coupling limit to arrive at dynamical 
semi-groups [29 1). In our case this is necessary because 
the field operators contain a factor of y/e. In the continuum 
limit where n — )• oo and e — )• 0, we arrive at a differential 
equation for p{x) with x G [0, L], 



dp 
dx 



-i[K,p\ + ]^{{^^\R,[R,p]]+Kc) (3) 

; ((^^^)i^[i^^p] + (^^^)[p,i?^]i? + h.c.) , 



[K, p]-\iZ ([M]M,,p\^ - 2M,pM] 



where in this expression all operators are evaluated at po- 
sition X, i.e., K = K{x), ^ = ^(x), etc. This is an 
example of a generator of dissipative dynamics which are 
manifestly completely positive, and up to dependence of x 
being of the form of a Lindblad generator 1,29 J . Here, the 
Lindblad operators are identified as Mi = iaR — hR\ 
Ma = iaR + hR\ M3 = cR\ and M4 = dR, where 
q2 ^ {{^,^f)/2, 52 = ((^f)/2, c^ = (^t^)^ and d^ = 
^^^t^ We write this equation as p'{x) = Cx{p{x)). 

We now describe a key feature of cMPS, shared by MPS 
and PEPS |30|, namely, their holographic property. What 
we mean here is that is the dynamics of a one-dimensional 
quantum field theory ^ in a cMPS is described by the dis- 
sipative dynamics of the boundary zero-dimensional field 



theory B alone. This is strongly reminiscent of the holo- 
graphic principle 11311 . The first step towards the holo- 
graphic property is to show how expectation values of field 
operators may be obtained in terms of the dynamics of 
the auxiliary system B alone. This is easy to establish 
using the following calculational principles, introduced in 
IT31 . Let A be any observable on A which is some prod- 
uct of the field operators and their derivatives at locations 
Xi, . . . ,x„ £ [0, L] in the continuum. The first step is 
to put the observable into normal order with all field an- 
nihilation operators on the right. Now we must calculate 
tv[AaA] = tr^ [{A(^Ib)U{L){uj_a0p{O))W{L)]. To 
eliminate the field operators we exploit the formula 

['^{x),U{L)] = -i f U{L-s)['^{x),F{s)]U{s)ds, 



where F{s) = Q{s) + R{s) ^ ^^(s), to commute 
all field annihilation operators past U{L). Temporarily 
assuming f3^ that ^(x)\i})^ = we thus learn that 
'^{x)U{L)\n)^ (g)I = U{L - x)R{x)U{x)\n)_A 
I. Similarly, to evaluate derivatives of ^ we follow the 
same procedure with an additional integration by parts to 
eliminate the derivative of the delta function. We find 

^'{x)u{L)\n)^mB = £{u{L-x)R{x)u{x))\n)_A(^ 

Ib = U{L-x){-[Q{x),R{x)] + R'{x))U{x)\n)_AmB- 
Higher derivatives may be evaluating using these methods, 
but are considerably more tedious. To proceed, it is expe- 
dient to employ the Jamiolkowksi isomorphism where op- 
erators M = J2ik "^jfcli)(^l ^^^ identified with quantum 
state vectors via \M) = J2ik ^-fclji k) ||35l . E.g., in this 

way, Eq. ^, with (^^^) = {^^) = (^^^) = 0, and 
(xl/ij/t) = lis written as 

^\p{x))=L{xMx)), (4) 

where the L/oMv///jfln is now L(x) = —iK0l+iI(^K'^ — 
|(i?^i?(g)I-2i?(^i? + I(g) R'^R) . Thus, to evaluate a 
correlation function tr[Ao"^] we simply have to integrate 
Eq. (Wji with additional insertions of the operators i?(x)(^I 
at the locations of ^, 10 R{x) at the locations of ^\ 
{-[Q{x),R{x)] + R'ix)) (g) I at the locations of ^', etc. 
We have thus completely eliminated the field A. 

To complete the derivation of the holographic property 
we need to show how to differentiate a one-parameter fam- 
ily of a cMPS (T^(t), where t is (real or imaginary) time. 
This is easily achieved using -j^U{L,t) = /^ U{L — 
s){§iF{s, t))U{s)ds, which implies that ti{Aa'^{t)) can 
also be evaluated in terms (an integral) of the solution of 
the Lindblad equation with an insertion of {-^F{s,t)) (g) 
I + I(g)(fF^(s,t))ats = t. 

There is a particularly convenient way to package the 
procedure we have just described: we define the generating 
functional Z[J], 



Z[J] = ([7|rexp 



dx Lx + J(A(x), /u(x)) 



Uo 



m?) 



where J{X{x), n{x)) = X{x)R{x) 1 + A(x)I (g) 
R{x) + fj.{x){-[Q{x), R{x)] + R'{x)) ® I + Jl{x)I (g) 
{—[Q{x), R{x)] + R (x)). Using Z[A(x),^(x)] we can 
obtain the expectation value of any field operator ^ or 
its derivative ^' via functional derivatives with respect to 
A(x) and ^{x) (see, e.g., chapter 10 of Ref. [36i for a brief 
introduction to functional derivatives), e.g., 



[■^\z)-^'{y)-^{x)) 



53 Z [J] 



5X{x)5^l{y)5\{z) 



(6) 



A,/j=0 



Now we have introduced the generating functional Z[J] it 
is straightforward to generalise the cMPS construction to 
arbitrary zero-dimensional field theories (i.e., D — )• 00) 
via the path-integral prescription: we promote K and R to 
functions of field operators and evaluate the time-ordered 
integral Eq. (j5]l via the standard path-integral prescription. 
Since the boundary field theory is a non-equilibrium theory 
it is convenient to employ the Keldysh formalism to eval- 
uate the generating functional. In this way we can easily 
understand some general properties of cMPS. 

For finite D, and in the translationally invariant set- 
ting where R and K do not depend on x, it is partic- 
ularly transparent to see that generically, all correlation 
functions decay exponentially. Let us assume that the Li- 
ouvillian L generating Markovian dynamics has a unique 
zero eigenvalue and the real part of any other eigenvalue 
is bounded from above by —A, A > being a gap. 
Since for ^(x)|i7)^ = we have that (^■f(xi)^(x2)) = 
iv^[U{L-X2)RU{x2){uJA®PB)U^{x^)R^U\L-x^)l 
we can apply the above rule for integrating the master 
equation, using techniques to compute two-point correla- 
tion functions for dynamical semi-groups |29|, to see that 
there exists a suitable c > with |(*I'^(xi)^(x2))| < 
gg-A(a:2-xi) Similarly, all other spatial correlators of our 
original field A are clustering. 

The perspective offered here from the viewpoint of con- 
tinuous measurement and the holographic principle allows 
us to easily generalise the cMPS ansatz class to field the- 
ories with a larger geometric dimension. We describe this 
generalisation in terms of a 2 + 1-dimensional bosonic the- 
ory, because the generalisation to higher dimensions of- 
fers no new complications (such a generalisation was an- 
ticipated in [,15J ). Intuitively our construction may be de- 
scribed as follows. Suppose we have a two (spatial) dimen- 
sional bosonic field theory A with field operator ^(x, y), 
x,y G [0,L]. Motivated by our 1 + 1-dimensional con- 
struction we propose to introduce an auxiliary 1 + 1- 
dimensional field theory B described by a tuple of field 
operators <!>„(?/), a = 1,2, ... ,D, (which may be trans- 
forming as a spinor, vector, or matrix etc.), which we think 
of as living vertically "at the boundary" of A. To prepare 
a quantum state for A we work entirely analogously as to 
before: we initialise A in some fiducial state, say the vac- 
uum I $7). We then interact an "infinitesimally thin" vertical 
strip of A and B according to some spatially local interac- 




FIG. 1. Here we illustrate the physical process underlying the 
construction of a 2 + 1-dimensional generalisation. The system A 
is initialised in the vacuum state jfi) and then infinitesimally thin 
vertical strips at horizontal location x are sequentially interacted 
with B at time t = x. 



tion 7^(0), where 

n{x) = i I dyR^{<^^{y),<^'^{y))(^^\x,y)+h.c., 
.Jo 

(7) 
and where Rx{^aiy),^'aiy))^ which may depend on 
the position x, is some polynomial in the field operators 
^a{y), their derivatives, and their adjoints at location y. 
We then proceed by interacting infinitesimal vertical strips 
at locations x of ^ and B sequentially at times t = x. In- 
terleaved between each interaction between the strip at x 
and B we also evolve B according to free dynamics gener- 
ated by some 



/C(x) 



dyK,{<^^iy),'l>'^iy)), 



(8) 



where K^ (<I>q (y) , <I>^ (y ) ) is some hermitian operator writ- 
ten in terms of a polynomial in the field operators ^^iy), 
their derivatives, and their adjoints at location y. This uni- 
tary process, illustrated in Fig. (fl), is described by the prop- 
agator 



U(L L) = ■j-Q-\fo dx ic(x)+{m{x)+h.c.) ^ 



(9) 



We can now apply, without modification, the analysis de- 
rived earlier to describe the holographic property of the 
states generated by U{L, L). A central role is again played 
by the Lindblad equation 



dp 
dx 



(This is the equation for the case where the initial state of 
A is the vacuum. The case of finite filling is more involved 
but is of a similar form to the 1 + 1-dimensional case.) This 



-i[IC{x),p] 

dy [Rl{^^{y)),R,{^^{y))p] + h.c. (10) 



describes a local dissipative field theory for B; expectation 
values of physical operators may be recovered by integrat- 
ing this equation with the appropriate insertions. Note that 
since the entire process is unitary all expectation values 
arising from integration of this equation are physical. 

While it is clear that, in the limit L — )• oo, when K^ 
and Rx are translation invariant the states arising from this 
construction will be translation invariant [i33;l it is far less 
obvious what conditions must be imposed on K^ and R^ 
in order that the resulting state is rotation invariant. This is 
a more subtle problem and there are two points of view one 
can take here: the first is that we simply ignore the rotation 
invariance issue and when we use 2D generalisation as a 
variational class we assume that the optimal state will in- 
herit the symmetries of the hamiltonian. The second point 
of view is that we should explicitly characterise the rotation 
invariant instances and use the resulting subclass as a vari- 
ational class for rotation-invariant systems, as this ought to 
be a simpler variational problem. In both cases we need to 
at least satisfy ourselves that rotation invariant 2D gener- 
alisations exist. In order to argue this we suppose that the 
dissipative boundary theory for B has a unique spacetime 
rotation invariant fixed point, which is equivalent to trans- 
lation invariance and the condition that Rxi^aiu)) trans- 
forms as a scalar under spacetime rotations. Then all we 
need is that the free dynamics of B generated by /C is also 
spacetime rotation invariant. This is a nontrivial condition 
and requires that the field ^^ transforms in a representation 
of SO{2), i.e., as a vector. These two conditions suffice for 
the resulting 2D generalisation to be rotation invariant. In 
the same way: if we want that the 2D generalisation is con- 
formally invariant then we require that the boundary theory 
for B is also conformally invariant ll37l . 

In this work we discussed an interpretation of a recently 
introduced variational class, continuous matrix product 
states, for bosonic quantum fields. We have explained how 
this class arises naturally from the procedure of continu- 
ous measurement, and used this observation to explain the 
key physical properties of cMPS, including the clustering 
of correlations. We also discovered a fundamental holo- 
graphic property possessed by cMPS, namely that the dy- 
namics of a quantum field in an cMPS can be completely 
understood in terms of a (dissipative) boundary field the- 
ory of one dimension lower. Finally, we have expressed 
the definition of an cMPS in purely field-theoretic terms 
which facilitates generalisations to higher-dimensional sys- 
tems and the use of perturbation theory. 

We have also pointed out that the zero-dimensional 
boundary field corresponds to the internal degrees of free- 
dom of an atom in cavity QED experiments; this opens up 
the possibility of simulating quantum field theories with 
simple dissipative dynamics. 

A great many future directions present themselves at this 
point: one can explore the utility of cMPS as a variational 
ansatz for numerical calculations. Additionally, the exten- 
sion of the theory of cMPS for fermionic fields and gauge 



fields is possible llJTl . Finally, the relationship of the reg- 
ulator that cMPS provides with other standard renormal- 
isation prescriptions (e.g., dimensional regularisation) re- 
mains to be elucidated. 
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state, which is a coherent state. Supposing the system is in 
1 0) the effect of measuring R can be effected by coupling 



the system and meter according to the interaction 



t_,-c.t, 



Ht =iR®a' -iW ®a 



(11) 



and evolving for a time t. If we temporarily assume that R 
is a normal matrix, [_R, R*] =0, then the resulting state vec- 
tor of the system and meter is given by e~'*^|*I')|0) = 

Ejii'/'jki)kj*). where R = J2f=irj\rj){rj\ is the 
eigenvalue decomposition of 7?, r-j are the complex eigen- 
values of R, and \a) — D{a,a)\Q) is a coherent state 
vector centred on the phase space point a G C. If we 
increase i we can arrange for \rjt) to become arbitrarily 
well distinguished by a physical measurement of the me- 
ter alone (eight-port homodyne detection |,34J ). Note that 
this procedure is physical even though the measurement of 
a non-Hermitian operator does defy immediate physical in- 
terpretation: it is a result of unitary evolution and standard 
quantum-mechanical measurement. 
[28] This is related to an interpretation supplied in ifTSl . 
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